Inflationary preheating and primordial black holes 
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Preheating after inflation may over-produce primordial black holes (PBH's) in many regions of 
parameter space. As an example we study two-field models with a massless self-interacting inflaton, 
taking into account second order field and metric backreaction effects as spatial averages. We 
find that a complex quilt of parameter regions above the Gaussian PBH over-production threshold 
emerges due to the enhancement of curvature perturbations on all scales. It should be possible to 
constrain realistic models of inflation through PBH over-production although many issues, such as 
rescattering and non-Gaussianity, remain unsolved or unexplored. 
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Introduction - The issue of whether initial conditions 
at the Planck era were suitable for the onset of infla- 
tion is both complex and controversial With these 
subtleties aside, there remains a cavernous space of pos- 
sible inflationary models |^ . The requirement of graceful 
exit from the cold inflationary phase into an acceptable 
radiation-dominated FLRW universe has proven a pow- 
erful filter on this model space. 

Failure to exit gracefully spelt the end of the old infla- 
tionary scenario H, is perhaps the major stumbling block 
in pre-big-bang models |^ and continues to plague string 
and supergravity models of inflation through the threat 
of overproduction of dangerous relics such as moduli and 
gravitinos Q. 

Perhaps the most radical way to end inflation is via 
preheating (see e.g., iQ]) in which runaway particle pro- 
duction occurs in fields coupled non-gravitationally to 
the inflaton. This explosive growth of quantum fluctu- 
ations drives similar resonances in metric perturbations 
on scales which range from cosmological to sub-Hubble 

i- . 

It is now recognized that in certain models preheating 
can alter the predictions of inflation for the Cosmic Mi- 
crowave Background (CMB) [^-13| by exponentially am- 
plifying super-Hubble metric perturbations. This does 
not violate causality but depends sensitively on the pre- 
ceding inflationary phase which determines the spectrum 
of X fluctuations |14|~fi8| . In this paper we discuss what 
appears to be a more robust mechanism for constrain- 
ing models of preheating - over-production of primordial 
black holes (PBH's). 

The idea that the amplification of metric perturba- 
tions during preheating would lead to enhancement of 
PBH abundances was raised early on and has been 
alluded to frequently since; e.g., Recently Green 

and Malik have used a semi-analytic approach which 
incorporates second order x field fluctuations to study 



PBH formation in a two-field massive inflaton model. 

Their results suggest that during strong preheating 
(<? ^ 1 0): PI3H formation could violate astrophysi- 
cal limits before backreaction ends the resonant growth 
of X fiuctuations. This is a crucial issue since strong pre- 
heating is generic in many models of inflation. However, 
Green and Malik used the results of 0| for the estimate 
of the time at which backreaction ends the initial reso- 
nance. As they point out this estimate does not include 
metric perturbations or rescattering and hence could be 
misleading. 
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FIG. 1. A schematic figure showing the numerical ap- 
proaches to preheating with numbers in brackets denoting 
appropriate references. See the text for discussion. 

Here we present first estimates of PBH production in- 
cluding backreaction computed dynamically. We find 
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that while preheating may lead to over-production of 
PBH's in some regions of parameter space, the result is 
sensitive to many subtle issues. 

To place our methods in context, consider Fig. 1 which 
shows the different numerical studies of preheating un- 
dertaken in the literature. The eventual goal of these 
studies is a fully nonlinear analysis of multi-field pre- 
heating including metric perturbations. So far this has 
been achieved without metric perturbations ("no $") - 
often with simplified expansion dynamics - through lat- 
tice simulations pi] ]. The furthest the community has 
progressed | |l9[ | in solving the full Einstein field equations 
is in a model with plane wave symmetry and a single 
scalar field. 

An alternative to full lattice simulations of preheating 
is the use of the Hartree, large-N, and mean field approx- 
imations ||2^ . Recently the Hartree approximation has 
been combined with the linear approximation for metric 
perturbations $ [^2|jl|JT|l and, in with the sec- 

ond order metric perturbations formalism of Abramo et 
al. It is this latter approach that we adopt. 

Immediate goals are fully nonlinear spherically sym- 
metric simulations suitable for studying individual PBH 
formation (c.f. j24j) and inclusion of rescattering effects 
in the presence of metric perturbations, <f>. The latter re- 
quires going beyond the Hartree approximation and eval- 
uating double and triple convolutions. 

The Model - We consider the two scalar field chaotic 
inflation model 



calculate PBH constraints in terms of the curvature per- 
turbation C.k rather than ^k- Ck is defined in terms of <I>fc 
and the Hubble parameter, H = a/ a, by 



(0.3) 



and is usually conserved on super-Hubble scales in the 
adiabatic single field inflationary scenario. In the multi- 
field case which we consider in this paper, this quantity 
can change nonadiabatically due to the amplification of 
isocurvature (entropy) perturbations. 
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where (j) is an inflaton field. During inflation x decreases 
rapidly towards zero if /\ ^ 1 in which case the 
temperature anisotropics in the CMB simply scale as 
AT/T ^ \/X. We therefore choose a self-coupling of 
A = 10^^^. During preheating, x ^-nd 5xk grow expo- 
nentially in very specific geometric channels or resonance 
bands which are well understood in terms of Floquet the- 
ory 110. 

We assume a flat background FLRW geometry with 
perturbations in the longitudinal gauge |8| : 



(1 + 2^)dt^ + a^{l ~ 2^)5ijdx'dx\ 



(0.2) 



where $ = "I>(x, t), the natural generalization of the New- 
tonian potential, describes scalar perturbations and a — 
a(t) is the scale factor. We decompose the scalar fields 
into homogeneous parts and fluctuations as 0(i, x) — > 
(t){t) + Scl){t, x) and x{t, x) ^ x{t) + ^x{t, x). 

The structure of the linearized Einstein field equations 
for this system can be schematically written in terms of 
two vectors: one for the FLRW background dynamics 
X = (0, 0, X, X, a, a), and one for the perturbation vari- 
ables in Fourier space: = {5(j)k,S^k,Sxk,SXk,^k,Ck)- 

While we solve the system of linearized Einstein field 
equations in the longitudinal gauge, it is convenient to 
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FIG. 2. An illustration of primordial black hole (PBH) for- 
mation during preheating due to growth of density perturba- 
tions. The PBH event horizon is schematically shown by the 
white ring in the final panel. Astrophysical limits on PBH's 
constrain /3, the ratio of PBH to total energy density. To con- 
strain theory one needs to map /3 into the mass variance a, 
most easily achieved with a Gaussian or chi-squared assump- 
tion for density perturbations. It is a that we calculate in our 
simulations. 

We include backreaction effects to second order in both 
field and metric perturbations pst , which implies that 
we integrate coupled integro-differential equations. The 
precise structure of these equations and additional de- 
tails can be found in Appendix and |]l^- [r^ , p3| . Here we 
illustrate the skeletal structure of the system, which has 
the form 

X = F(X,(Y2)), 
Yfc = G(X, {Y'))Yk, 
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pert; 
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where the variance is defined by 
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(0.5) 



for any field o. F and G are nonlinear functions of 
the spatially homogeneous background vector X and the 
variances of the components of Y. The complete system 
is integrated from 50 e-folds before the end of inflation 
to provide the appropriate initial conditions for preheat- 
ing. The initial values at the start of inflation are chosen 
as ^ = 4mpi and x — 10^'^mpi with conformal vacuum 
states for the fluctuation^. Including the field variances 
ensures total energy conservation at 1-loop. 

Primordial black hole constraints - Since PBH's form 
from large density fluctuations |2^, it is an obvious 
concern that preheating might encounter problems with 
PBH constraints arising from the Hawking evaporation of 
small PBH's or from overclosure of the universe (ripBH > 
1) for heavy PBH's. 

To quantify this suspicion one needs to compute the 
mass function /? (I^S: 



Ppbh 
Ptot 



P{S) dS, 



(0.6) 



where P{5) is the probability distribution of the density 
contrast, (5, and dc, (~ 0.7) is the critical value at 
which PBH formation occurs in the radiation dominated 
era. 

Usually one assumes a Gaussian distribution P{S) — 
l/(\/27rcr) exp[— (5^/(2(T^)], where a is the mass variance 
at horizon crossing. Observational constraints imply that 
P < 10~^° over a very wide range of mass scales, which 
translates into a bound on the mass variance of cr < cr* = 
0.08. (J > cr* corresponds to PBH over-production in 
the Gaussian distributed case. When the distribution 
is instead first order chi-squared - an approximation to 
the X density fiuctuations in preheating (see the later 
discussions) - the threshold is = 0.03 ||2(]|| . 

Defining the power-spectrum of the curvature pertur- 
bation as V( = fc'^|CfcP/(27r^), the mass variance can be 



expressed as ^0 
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A) r,wikR)f. (0.7) 



We choose a Gaussian-filtered window function W{kR) = 
exp(— fc^i?^/2) where R = is the artificial smoothing 
scale [Q. We can expect exponential increase of a due 
to the excitement of field and metric perturbations dur- 
ing preheating. We solved the Einstein equations (0.4) 
numerically, varying the ratio and evaluated the 

mass variance with two cut-offs fc, = aH and fc» — lOaH 
to investigate sensitivity to cut-off effects. 
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FIG. 3. Threshold PBH formation - the growth of a, 
C,k = k^^'^Ck, and = k^^^&Xk/nipi for a super-Hubble 
mode K = k/{\/\(f>o) = 1Q~^^ vs dimensionless time 
X = \/\<f)orj in the case g^/\ = 2.5, where <^o — O.lrripi 
is the value of inflaton when it begins to oscillate coher- 
ently. With the choice fc, = aH in the window func- 
tion W{kR) — exp{—k^R^/2), a just reaches the threshold 
CT, = 0.03 for the PBH formation for chi-squared first order 
distributions. 

When Xk fluctuations are amplified during preheating, 
this stimulates the growth of the metric perturbation, ■ 
On cosmological scales this effect is sensitive to the sup- 
pression of X and Sxk modes in the preceding inflationary 
phase. 

When g^/\ = 0(1), this suppression is weak since the 
X field is light and once the long-wave Sxk modes 
grow to of order S(j)k during preheating, super-Hubble 
and Cfe are amplified until backreaction effects shut 
off the resonance. This amplification occurs in the region 
1 < g'^/X < 3 where the A; ~ modes lie in a 

resonance band. The increase in Cfc leads to a correspond- 
ing growth of the mass variance a which can reach the 
threshold cr* = 0.03 for 1 < g^/\ < 3 and 6 < g^/X < 10 
with the cut-off set at fc* ~ aH, i.e., around the Hubble 
scale (see Fig. ^) . 

As is increased, the x field becomes heavy and 

suppressed during inflation. This restricts the amplifi- 
cation of super-Hubble metric perturbations [ p^ despite 
the fact that the fc ^ mode of Sxk lies in a resonance 
band for n(2n-l) < gVX < n{2n + l), n = 1,2,3... [|5), 
as is evident from Fig. |j. However, since sub-Hubble Sxk 
modes are not suppressed during infiation 1 14 |] , and 
C/c on sub-Hubble scales do exhibit nonadiabatic, reso- 
nant, growth for g^/X ^10, which leads to growth of 
a. 

However, we do not find that this is significant enough 
to lead to cr > cr* for g^/X S> 1, except for very short 



*Using the initial condition x = 10 ^mpi we reproduced the 
results of Ref. [||. 



'We have reproduced the result that the homogeneous part 
of the X field is ampUfied by the second order couplings be- 
tween <l>fe and Sxk [M despite of the infiationary suppression. 
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intervals around = 0, in contrast to the expectations of 
[^o| . However, when we enlarge the cut-off frequency 
fc* to IQaH, we do find cr > 0.08 in wide ranges of 
parameter space (see Fig. |^). Somewhat surprisingly, 
these super-threshold regions arc all clustered around the 
super-Hubble resonance bands in /X space. 
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FIG. 4. The power spectrum of C, at the end of preheat- 
ing for the values of g^/X = 2,50 with fc« = aH. Since the 
inflationary suppression of x becomes stronger as /\ is in- 
creased, the growth of C,k at long wavelengths is suppressed. 
Note also the dominance of the (?^/A — 2 modes at sub-Hubble 
scales K > 1. Inset : The evolution of 5xk for a super-Hubble 
mode K = k/{VX(j)o) = 10"^^ for g^/\ = 2,50. The suppres- 
sion of the initial conditions, due to the preceding inflationary 
phase, in the heavy case g^ /X = 50 is evident. 




FIG. 5. The mass variance a vs g^ /\ for two window func- 
tion cut-offs fc* — aH and fc* = 10a//. The threshold of 
(7, — 0.03 in the chi-squared distributed case is shown and 
is marginally crossed for the regions around g^/X ~ 2 and 
g^/\ ~ 8 when fc* = aH. For fc* = 10a// a quilt of regions 
above the Gaussian threshold cr, = 0.08 emerges which coin- 
cide closely with g^ /X — 2t? corresponding to Floquet indices 
with maxima at longest wavelengths. 



In our model these fields are represented by x ^-^d fhe 
initial value is set 55 e-folds before the end of inflation. 
This problem has two facets - the initial value of the 
background, or vacuum expectation value of x, and the 
initial value of the distribution of fluctuations, i.e., 5xk- 

A sensible choice for the latter is the Bunch-Davis vac- 
uum, but it is the initial value of the homogeneous part 
of X, which is of the most importance, since if x = 
(the minimum of the potential) no resonance can occur 
at linear order. 

We have found four suggestions for setting the initial 
value, Xii follows. 

(1) Choose the value of x which maximises the proba- 
bility distribution in eternal inflation for fixed large val- 
ues of the inflaton (0 > Impi) at a specific time. Since the 
regions with the largest Hubble constant dominate the 
distribution this corresponds to choosing x > l"^pij 
i.e., super-Planckian chaotic initial conditions for x- 

This suggestion is, however, sensitive to the choice of 
a hypersurface for setting the initial conditions. If one 
defines initial conditions on the hypersurface of energy 
density equal to the Planck energy for instance, then the 
Hubble constant will likely be maximised by placing all 
energy into the field with the flatest potential, rather 
than distributing it amoung the various fields, some of 
which may have steep potentials. This will lead to van- 
ishingiy small initial x unless x is a good infiaton, i.e., 
.9' < ^• 

(2) Choose Xi to satisfy xf = (x^) ||l 
Bunch-Davis condition, Xk ^ '^1 ^ 
estimate (x^) = 2^ / k^\Xk^dk to be 
variance is super-Hubble dominated during inflation and 
where k^, = H is the natural cut-off at the Hubble scale. 

Now if inflation is driven by then ~ V{(j))/mp^ 



|. If we use the 

l/^/gcj), we can 
^ k^Jig^) if the 



and we find xf 
(I 



l/(0)^/^/(50TOpi). For the potential 
1|) this leads to the estimates Xi ^ 10~^mpi for g^/A = 
0(1) and Xi ~ 10"*g"^/^TOpi for g^/X > 1 if we take 
A 10"^^ and (j)^ ^ 4mpi. 

(3) Choose the value of Xi which leads to a station- 
ary distribution in eternal infiation (where the classical 
drift and quantum fluctuations are balanced) ^. Assum- 
ing quantum fluctuations S(j) ~ H/2t: on characteristic 
time scales 6t ^ H^^ one arrives at Xi ^ /{g^^"^) ^ 
hence y, ^ 10"'' 



and 



- IQ- 



'pi 



for 



A3/^<^V(5^-^i) 
gVA = 0(1) and x^ ~ 10"^^5"^mpi for g'^/X > 1. 

(4) Finally we may choose the value of Xi which cor- 
responds to the instantaneous minimum of the potential. 
It suggests Xi = 0- This argument has several problems 
the most fundamental of which is that the system is not 
in equilibrium since the x field is not strongly coupled 
except for g^ /X 3> 1. 

Despite the wide range of possible initial values, Xii 



Initial conditions for the x field - An issue of general 
importance which has been little studied is that of initial 
conditions for non-infiaton fields at the start of inflation. 



*We thank Alan Guth for this suggestion. 
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the start of inflation, Fig. || shows that the final mass vari- 
ance, and hence the probability of PBH over-production, 
depends rather weakly on Xi- 



0.010 
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FIG. 6. The smoothed dependence of the final mass vari- 
ance a on the initial condition, Xi, 55 e- folds before the start 
of preheating for g^/X — 2,50. Note the relatively weak de- 
pendence in both cases. 

Potential problems and unresolved issues - Our results 
suggest that PBH over-production may not be generic in 
strong preheating. However they can only be considered 
as preliminary for a number of reasons: 

• There are at least two fields critically involved in 
preheating. Even if the inflationary fluctuations 
are Gaussian, the fluctuations induced by preheat- 
ing are typically not. If the x held has no vac- 
uum expectation value, its density fluctuations are 
roughly oc 5x^ , so approximately chi-squared if Sx 
is Gaussian distributed. As discussed above, we 
take X 7^ 0. The recent results of suggest that 
X is Gaussian distributed before rescattering sets 
in and hence the density perturbations would be 
Gaussian, at least while dominated by linear fluc- 
tuations. 

Rescattering leads to non-Gaussian distributions 
and to S(j) oc Sx'^ 0- The applicability of the cri- 
terion cr > (7* therefore depends largely on when 
PBH formation actually takes place - before or af- 
ter rescattering. 

Further, the density fluctuations may go nonlinear. 
Since 6 e [— l,oo) this necessarily skews the distri- 
bution, similar to the toy-model discussed in the 
second Ref. of j2^. Non-Gaussianity may drasti- 
cally alter the relationship between (3 and a |^,^ , 
changing cr* and requiring the use of higher-order 
statistics. 

• In preheating the Hubble radius is vastly smaller 
than the true particle horizon and resonance bands 
often cover the complete range of scales. Predicting 
the mass spectrum of PBH's created during pre- 
heating is therefore a subtle issue. Grudely one ex- 



pects a wide range of PBH masses to be produced, 
even without criticality arguments |^ . 

This is related to our results showing cut-off, k^, 
sensitivity. The increase in a when is altered 
from aH to lOaH reflects the important contribu- 
tions of sub- Hubble modes. Does this necessarily 
imply that the resulting PBH's are very small ? If 
so, they are not constrained since they evaporate 
harmlessly long before nucleosynthesis. 

• We have not included rescattering. This is known 
to enhance variances over the Hartree approxima- 
tion at small resonance parameters, q, in the ab- 
sence of metric perturbations For q ^ 1 
however, the situation is reversed and variances 
are overestimated by the Hartree approximation. 
Whether these results are stable to inclusion of met- 
ric perturbations is unknown, but this may provide 
a way to avoid PBH over-production since it should 
filter through to Cfe and a. 

• Fig. H shows fj as a function of g'^/X. The value of 
a plotted is its maximum at the end of preheating. 
However, a does grow larger than this value, instan- 
taneously exceeding O.OI, even for fc* = aH, when 
4> — 0. We choose the more conservative route of 
not taking these as the true maxima, but the ques- 
tion remains, can large a, attained for very short 
periods, lead to PBH formation ? 

• We solved the x field equation, including second 
order terms such as ($(5x) fl^- Initially Sx and 
$ are correlated but when the x fluctuations are 
sufficiently amplified, they are well described by 
classical stochastic waves which may be un- 
correlated with metric variables. 

It is uncertain that contributions of second order 
metric terms should be included during such classi- 
cal regimes. Since this issue affects x rather signifi- 
cantly, the quantum to classical transition appears 
to be of quantitative importance, deserving further 
study. 

Conclusions - We have studied primordial black hole 
(PBH) formation during preheating using numerical sim- 
ulations of the perturbed Einstein field equations includ- 
ing second order field and metric backreaction effects. We 
found that there exist parameter ranges where standard 
Gaussian and chi-squared thresholds for PBH formation 
are exceeded. 

Nevertheless, the results are not unambiguous. We 
discovered a significant sensitivity to the window function 
cut-off, fc*, and since preheating is expected to lead to 
non-Gaussian fluctuations, it is not clear how realistic the 
Gaussian threshold for PBH formation is. Nevertheless, 
PBH over-production constraints are very robust. The 
study of PBH's in preheating is an exciting area which 
may lead to strong constraints on realistic inflationary 
models. 
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We note that there are a number of possible es- 
cape routes to preserve preheating but avoid PBH over- 
production. Fermionic preheating is very unhkely to lead 
to PBH formation unless the fermions are extremely mas- 
sive. Similarly, instant preheating p^ , which draws en- 
ergy away from the x held almost immediately seems 
likely to stall PBH formation, as does a large x self- 
interaction. 

On the other hand, since growth of (k and a is seeded 
through isocurvature/entropy perturbations [l^ , it is 
possible that other models of reheating, such as non- 
oscillatory models , which lead to significant isocurva- 
ture modes, may also have a PBH over-production prob- 
lem. 

Nevertheless, the precise scenario of the PBH forma- 
tion during preheating can only be understood properly 
by overcoming two serious hurdles - (1) understanding 
the probability distribution of density fluctuations during 
preheating and (ii) going to fully nonlinear simulations of 
resonant PBH formation which include rescattering and 
nonlinear metric perturbations. 

Acknowledgements - We thank Christopher Gordon, 
Alan Guth, Anne Green, Roy Maartens, Kei-ichi Maeda, 
Karim Malik, Masaaki Morita, and Takashi Torii for use- 
ful discussions and comments. 



2(^Sx) - 4($5x) - GH{<^>5x) 
2 



4x($$) 



(<i>v^(Jx)> = 0, 



(0.10) 



where G = m^^ is Newton's gravitational constant. Note 
that (...) implies a spatial average. In spite of the expo- 
nential suppression during inflation, operative when the 
X field is heavy (g^/A ^ 1), the x field can be signifi- 
cantly enhanced in the presence of th e second order met- 
ric backreaction terms in Eq. (3.1C) as pointed out in 
Ref. 0. 

The Fourier transformed, perturbed Einstein equations 



r,2 



-?,\{^'' + {5cp-'))+g\x^ + {5x^)) 
+ 2i4> + m.p)^k-2g^(l)xSxk, 



dxk + SHSxk 



g'iq^' + {5^')) 



5xk 



4x$fc + 2(x + 3i/x)*fc-2.g' 



(0.11) 



(0.12) 



APPENDIX: DETAILED FORM OF THE 
EVOLUTION EQUATIONS 

In this Appendix we present the evolution equations in 
details. We include second order field and metric back- 
reaction effects |p3| in the background equations, which 
is combined with the Hartree approximation p^ . 

Then the Hubble parameter and homogeneous parts of 
the scalar fields satisfy ||l^,0 : 



H' = 



inG 



1 :^ 



1 



1 



1 



{Sx' 



1 



{{W6xr 



2'"'^ ' ■ 2a2 
2(A03 + 52^x')($'50) + 2g'<i?x{^^x) 



2' 



+ 4i7($$) - ($2) + ((v$)2), 

3H^){1 + 4($2)) + A0((/)2 + 3((5(/)2)) 
.9'(x' + (<5x'))<^-2($(50)-4(<i>4) 



(0.8) 



6H{<^6<j>) + 4(/)($$) 



:($V2(50)) =0, (0.9) 



(X + 3Hx){l + 4(<i>2)) + ^2(^2 ^ ^s^2^)x 



$fc + = A7rG{(t>6(t)k + xSxk) 



(0.13) 



We find from Eq. ( 0.13| ) that metric perturbations grow 
if X and Sxk fluctuations are amplified during preheat- 
ing and the x-dependent source term exceeds the 0- 
dependent one. When field and metric fluctuations are 
sufficently amplifled, the coherent oscillations of the in- 
flaton condensate, </>, are destroyed. The entire spectrum 
of fluctuations typically moves out of the dominant res- 
onance band and the resonance is shut off. 
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